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Abstract: We perform lattice simulations of two-dimensional supersymmetric Yang-Mills 
theory with sixteen supercharges with a lattice action which has two exact supercharges 
(Sugino lattice action). According to the gauge/gravity duality, the theory at finite tem¬ 
perature is expected to be well described by the corresponding black 1-branes, at low 
temperature in the large N limit. We aim to confirm the duality conjecture by comparing 
the lattice results with the theoretical predictions obtained in the gravity side. In this 
article, at the beginning of this study, we examine the supersymmetric Ward-Takahashi 
identity to test whether the lattice action reproduces the correct continuum theory. Nu¬ 
merical results of the SUSY WTI strongly suggest us that any cut-off effects, which break 
supersymmetry, disappear in the continuum limit. In addition, we study the issue of de¬ 
generate vacua and find that the admissibility condition or any other constraints of the 
link fields which guarantee the unique vacuum are not always needed. 
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1 Introduction 

The correspondence between gauge theory and the theory of gravity is an important subject 
in theoretical physics. It was originally proposed as the AdS/CFT correspondence]!] which 
was predicted in string theory. Once we accept the duality, there are many interesting 
applications [2], The duality conjecture gives us deep insights into some long standing 
problems, such as the black hole information loss paradox[3]. Moreover, there are many 
attempts to solve strongly coupled theories, like QCD, analytically in the gravity side [4], 
[5]. To ensure the validity of the applications, it is important to give strong evidence of 
the duality conjecture. In this study, we aim to verify it by using lattice simulations of 
supersymmetric Yang-Mills theories. 

Lattice theory provides a framework to study strongly coupled gauge theories. In the 
context of the gauge/gravity duality, gauge theory is typically a strongly coupled one and 
thus is a good target for lattice theory. We can test the duality conjecture by comparing 
lattice results in the gauge theory side with theoretical predictions in the gravity side. In 
principle, unlike localization methods[6] and other analytical treatments]?] for BPS states, 
lattice gauge theory allows us to test it for any physical quantities and independently of N 
and ’t Hooft coupling constant A. Although lattice theory achieves great success for Yang- 
Mills theories, the nature of supersymmetry algebra makes it hard to put supersymmetric 
Yang-Mills theories on the lattice. Therefore, to test the gauge/gravity duality for SUSY 
theories using lattice methods has been a difficult issue for a long time. 
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In the last decade, some developments of lattice supersymmetry changed the situation 
completely. Several lattice formulations of SYM theories were proposed in [8], [9], [10], 
[11], [12], [13], [14], [15], [16], [17] and enabled us to perform direct tests of the duality. 
In one dimensional SYM with 16 supercharges, some evidences of the correspondence be¬ 
tween the gauge theory and the dual black hole[18] have been obtained from numerical 
simulations [19],[20]. In two-dimensional J\f = (8,8) SYM, lattice simulations have been al¬ 
ready performed and the phase structure of the theory was investigated [21], However, this 
subject deserves more studies because some physical quantities, like the internal energy of 
the black l-branes[22], have not been estimated yet. Therefore, we perform lattice simula¬ 
tions of J\f = (8, 8) SYM in two dimensions by employing the Sugino lattice actionfll]. 

In this paper, as the first step of our study, we examine the behavior of the Sugino 
lattice action for A f = (8, 8) supersymmetric Yang-Mills theory in two dimensions. By 
examining the SUSY Ward-Takahashi identity, we confirm that supersymmetry is restored 
in the continuum limit. This paper is organized as follows: We briefly review the Sugino 
lattice model and discuss degenerate vacua in section 2. Then, we see numerical results 
of the SUSY Ward-Takahashi identity in section 3. The last section 4 is devoted to a 
summary of our results. In appendix A, we give a relation between the notations used in 
this paper. The full lattice action can be found in appendix B. 

2 Two-dimensional J\f = (8, 8) SYM 

In this section, we explain the continuum action of the target theory and the lattice action 
we employed. We use Sugino’s method to define the lattice actionfll]. After introducing the 
continuum action we rewrite it as a Q-t-exact form, then introduce the lattice action from 
the Q-|--exact one. This way, the supersymmetries associated with the two supercharges 
Q± are kept exactly on the lattice. 

2.1 Continuum theory 

Two-dimensional supersymmetric Yang-Mills theory with sixteen supercharges contains 
two gauge fields = 0,1), and eight scalar fields Xi(i = 2, 3 • • • ,9), and sixteen fermions 
„(( y. = 1,2, • • • , 16) as the fundamental field variables. The scalars and fermions, as well 
as the gauge fields, are in adjoint representation of SU(N ) and can be represented as 
matrix-valued fields. 1 

The Euclidean action is given by 

S =jJ d ' X tT {\ F ^ + \ iD K X i) 2 ~ \\ X ^ X 3? 

+ U> a D^ a - JtfaMapDiVp + ^/ 3 ]| , ( 2 - 1 ) 

where, unless we use twisted fields which are defined later, the Roman indices i,j run from 
2 to 9, the Greek indices associated with the space-time directions /r, v take the value 0 

1 Let ip be a scalar field, a gauge field, or a fermi field. In this paper, ip can be expanded as ip = <p a T a 

where T a are group generators satisfying tr(T a T°) = 5 a b, 
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or 1, while the Greek indices associated with spinors a,/3 run from 1 to 16. The held 
strength is Fq\ = OqAi — 9\Aq + i[AQ,A{\ and the covariant derivatives are defined by 
Dpp = dptp + i[Ap, ip\. In addition, 7 a (a = l,--- ,9) are real symmetric matrices that 
satisfy the nine-dimensional Euclidean Clifford algebra, { 7 a, 7b} = 2 5 a b- 

This action is invariant under the following sixteen supersymmetry transformations, 


QaAp i(p)p)a/3^fh 

QotX{ — i(pyi)a/3^ 

Qo$ (3 = iMa P F 01 + ( 7 pXi + — (p/i'yj') a p[Xj i ,Xj\, 


( 2 . 2 ) 

(2.3) 

(2.4) 


where 70 = i. Furthermore, the two-dimensional Poincare group and the internal 0(8) 
rotations are symmetries of this theory. 

The easiest way to give the action (2.1) is to perform a dimensional reduction from ten¬ 
dimensional J\T = 1 SYM. However instead of going directly to two dimensions, reducing the 
theory to four dimensions first, we obtain four-dimensional A/" = 4 SYM which have a known 
topologically twisted version[23]. In the topological held theory, two of the supercharges 
Q± are scalars and nilpotent. Its action is written in a Q-(--exact form which is needed for 
the lattice action defined in the next section. Then, by reducing dimensions from four to 
two, we can give another representation of the action ( 2 . 1 ). 2 

As a result, the action ( 2 . 1 ) can be written in the Q-|--exact form with the four¬ 
dimensional topological held theory notation, 



(2.5) 


where /r runs from 0 to 3, and i,j, k from 0 to 2, and is a totally antisymmetric tensor 
that satisfy 6012 = 1- Note that p and i are different from those of the original variables, 
A^,Xi. Here, the bosonic helds consist of two gauge helds Ao,A\, six real scalar helds 
A2, A3, Bi,C , and two complex scalar helds (f>±. The fermions are given by ip±p, x±i- V±- 
The extended held strengths are 



( 2 . 6 ) 

(2.7) 


F\xv — dpA u d v Ap + 


with dp = 0 for /1 = 2, 3. 

The associated Q-(--transformation are 


2 


See appendix A for the relation between both notations. 
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( 2 . 8 ) 


Q±A^ — ip± ^ 

Q± B i Y : i ■ 
Q±C = rj±, 


Q±X±i = [ B u 0±], Q±X=r = ^[C) -Si] ± 
Q±V± = [CA±\, Q±V* = [0qF>^±], 


Q±</>± — 0) Q±<^=f — ^7=f ? 


Q±-Sj = ±([X=Fi,0±] + ^[x±i,C] + ^[-Si^i]): 

= ±([ip Tti , </>±] + ^[V’± M) ci - \ D ^±) ■ 


We introduced seven auxiliary fields H^ and H t to define Q± as closed transformations. 
With this arrangement, the Q-t-transformations are nilpotent up to gauge transformations, 
Q\ = id ( f >± and {Q + ,Q-} = — idc with 5^ being a gauge transformation with a gauge 
parameter ui, 


SojA/j, = -D^uj, 6 w <p = -i[tp,u;], (2.9) 

where p is either a scalar or a fermion. 

2.2 Lattice theory 

We consider the lattice theory on a two-dimensional lattice with a lattice space a, and 
hereafter we set a = 1 for simplicity. Let Nt,N x be the temporal and spatial lattice sizes, 
respectively. Then, the lattice sites x = (t, x) are labeled by integers, 

t = l,---,N t , x = l (2.10) 

The gauge fields are defined on the links as gauge group-valued fields (x ), and the other 
fields are defined on the sites. In this subsection, we impose periodic boundary conditions 
on every held for both temporal and spatial directions in order to discuss SUSY invariance 
of the lattice action. In section 3, the fermions satisfy anti-periodic boundary conditions 
for the temporal direction because we perform numerical simulations at finite temperature. 

The infinitesimal lattice gauge transformations for the link fields and the other adjoint 
fields ip are given by 

SujU^x) = iu](x)Un(x) - iUfj,(x)uj(x + fi), 6p(x) = - i[ip(x),u(x )\, (2.11) 

where fi is a unit vector in /i-direction and w is a gauge parameter defined on the sites. 
The forward and backward covariant difference operators are given by 

Vjv?(®) = U^{x)ip{x + p)Ul(x) - <p(x), (2.12) 

V~^(z) = <p(x) - t/J(x - $)ip(x - fi)Ufj,(x - p), (2.13) 


respectively. 

The use of the link helds leads to extra difficulties and forces us to adapt the lattice 
counterparts of the Qi-transformations in order to assure that they are nilpotent and to 
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keep gauge invariance intact. In fact, some modifications are needed only for the fields 
with directions n = 0,1. The modified (^-transformations for them are 


Q±i’±fi = -iV+<£± + ± M , 

Q±*h fv = ^v+C'±^ + |{V’+ m ,V’-m }> 

Q±H„ = ±( 4>± + + i |>± M , c + iv+c] 

- ^V+?]± + JtV’iM, V’-A*} ± . 


While the transformation law for the other fields is the same as the continuum one. From 
the definitions above, the (^-transformations are nilpotent, up to the lattice gauge trans¬ 
formations (2.11). We have Q\ = i5^ ± and {Q+,Q-} = — i8c even on the lattice[ll]. 

Let us give the lattice action from (2.5) by replacing the integral with a summation 
over the sites and changing the held strengths (2.7) to lattice counterparts, 


N 


t,X 


S = Q+Q-^^2 tr 1 ~ 4iBiF x - 3 e ijk BiBjB k 


1 


-V’+mV’-m ~ X+iX-i ~ 7’l+T f > 


where the extended held strengths must be carefully chosen, 


(2.15) 


F+ = I(V 0 + 7L 3 + V+7L 2 ), (2.16) 

F+ = i(Vr7l3-V 0 -yl 2 ), (2.17) 

F 23 = ^(*[^ 2 ,^ 3 ] +T 01 ), (2.18) 

with 

Foi(x) = - l - (p 01 ffl - Pji(x) - ^tr(P 01 (x) - 40?))) > ( 2 - 19 ) 

Poiix) = Uoffllhtf-+ 6)tf 0 (x + l)u\ (£). (2.20) 


See [11] for more details. If the forward-type difference operators are used for both (2.16) 
and (2.17), then the lattice action has fermionic and bosonic doublers. This mixed choice 
forbids the doublers. The hermiticity of Fq\ leads to a semi positive dehnite boson action. 
The trace-part of F 01 makes no contribution to the lattice action. Therefore we employed 
the traceless hermitian Fqi (2.19). Moreover, by performing the Q-(--transformation (2.14), 
a concrete form of the lattice action can be obtained. However, it is somewhat complex 
and lengthy, therefore we give it in appendix B. 

Note that, our simple choice of Toi (2.19) yields a problem of extra vacua, as was hrst 
pointed out in [9] and detailed in the next section. To avoid the problem, we can take 
an admissibility-type held tensor[10] or a tan(0/2)-type held tensor[13]. To use them is 
theoretically clear because the extra vacua are forbidden. But to write a simulation code 
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using such field tensors is complex and demanding task, therefore we chose the simple field 
tensor (2.19). As seen in section 3.3, no transitions from the desirable vacuum to the extra 
vacua are observed. Therefore, to use (2.19) is not a problem in practice, at least, for 
simulation parameters on which we focus in this paper. 

2.3 Continuum limit and degenerate vacua 

In this section, we discuss the naive continuum limit of the lattice action and the presence 
of degenerate vacua[9]. We argue that the degenerate vacua, while making the theoretical 
study of the continuum limit difficult, do not pose any serious problems in the simulations 
when testing the duality conjecture. 

Let us consider the naive continuum limit of the lattice action (2.15). The ’t Hooft 
coupling A has mass dimension two which gives the typical scale of this theory. Therefore, 
a 0 is realized by taking Ao(= Aa 2 ) —>• 0 with fixed A. Then, dimensions of the fields 
indicate 


U l “ L = 1 + iaA c ° nt ' + • • • , X^ at - = aX™ nt ', ¥ aL = a 3/2 ^ conL , (2.21) 

where the fields with lat. are dimensionless fields which define the lattice action (2.15), 3 
while the fields with cont. are the ones in the continuum action (2.1). After performing the 
(5±-transformations, the lattice action has many higher derivative terms. By neglecting 
the higher derivative terms and using (2.21), we find that, in the continuum limit, the 
lattice Q±-transformations coincide with the continuum ones and the lattice action (2.15) 
reproduces the correct continuum action (2.1). 

This argument of the continuum limit can be justified if 

P 0 i(x) -> 1, C 0, as A 0 -> 0, (2.22) 

where Poi are the plaquettes (2.20) and C is the sum of the higher derivative terms in 
the bosonic action. 4 When Pqi —1 1, the link fields also approach unity up to gauge 
transformations, therefore the expansion around one = 1 + iaA^ + • • • is justified. 5 
We only have to study the higher derivative terms C in the bosonic sector because, when 
C vanishes, the corresponding higher derivative terms in the fermionic sector also vanish 
thanks to the Q-t-symmetries. 

We have to examine (2.22) numerically for each simulation parameter because it does 
not hold in general as explained below. When the scalar fields are set to zero, the boson 

3 We can express the lattice action (2.15) in the original variables A )L . X L . by using the relations 
given in appendix A. 

4 The higher derivative terms in the bosonic sector are given by 

C = £ Y , tr {- X a] [W, X e ] + [X a , X+°] [X e , X+°] + [X + 1 , X e ] [X 3 , X+ 1 ] + ■ • • } , (2.23) 

t,X 

where X^ = Xi + VtXi. In the naive continuum limit, the three terms above cancel each other by the 
Jacobi identity. The complete equation of C is given in (B.16). 

5 Here, we assume that the lattice size is infinite. 
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action is given by the so-called double plaquette action, 


^lx<=o 


N 

2Aq 


E tr 




(2.24) 


We find that Ao —> 0 does not immediately mean Poi = 1- For example in the SU(2) 
case, P 0 i = — 1 is another possibility which gives a zero point of the integrand in (2.24). 
For non-zero scalar fields, the situation becomes more serious. For the simplicity of the 
explanation, let X 2 -,X% be non-zero and the other scalars be zero. Then, the boson action 
can be given as follows: 


where 


S B | 


A'i •—-X 9 0 


N 




(2.25) 


<A) = Vq"X3 + \ 7 fX 2 , (2.26) 

¥1 = Vj X 3 — V 0 X 2 , (2.27) 

(p2 = i[X 2 , -X 3 ] + Fbi- (2.28) 


Near the continuum limit, if V /x Xj ~ 0, that is, if Xj are smooth, then the higher derivative 
terms actually disappear. However, the conditions ipi = 0 imply that there can exist 
extra minima other than V^Xj = 0. For example, non-smooth configurations that satisfy 
V^Xs = —V+X 2 = 0(1) give (po = 0. Therefore, C could survive as Ao —> 0. 

Fortunately, the condition (2.22) is likely to hold for parameters in which we are 
interested. We actually consider the theory at finite temperature and test the gauge/gravity 
duality in the large N limit. We expect that the undesirable configurations are suppressed 
by temperature effects or by the influence of the large N limit. For the gauge fields, the 
two minima Poi = ±1 are separated by a finite barrier associated with the saddle point of 
(2.24). The height of the barrier becomes larger in the large N limit or in the continuum 
limit because the action has N /Ao as an overall factor. Furthermore, for the scalar fields, a 
potential reducing the fluctuations of the fields is dynamically generated by SUSY breaking 
temperature effects. At high temperature, the potential makes the scalar fields Xj stay 
near zero. These suppress transitions from the trivial vacuum = 1 and Xj = 0 to the 
others, for example, from Pqi = 1 to Poi = — 1 etc. Therefore, it is reasonable to consider 
that the lattice action with the naive field tensor (2.15) is sufficient to test the duality 
conjecture, despite the need of numerical confirmations of (2.22). 6 

6 As seen in section 3.3, starting the HMC from a trivial configuration 11^ = 1 and Xi — 0 at the 
temperature T/X 1 ^ 2 = 0.3 in SU(2), we observed that the higher derivative terms C were small in the run. 
Moreover, Poi tended to stay close to unity. In this case, we added a mass term for the scalar fields to the 
action. The mass term forbids large fluctuations of the scalar fields and contributes to the realization of 
C —^ 0 as Ao —I 0 for such low temperature and small N. While, somewhat surprisingly, Poi ~ — 1 was not 
observed for the simulation parameters even for aX 1 ^ 2 ~ 0.28. We have considered that the naive choice of 
the field tensor (2.19) works at high temperature, in the large N, or near the continuum limit. Contrary to 
our guess, it may work for a wider range of the simulation parameters. 
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3 Restoration of supersymmetry 

In this section we verify the restoration of the full supersymmetry in the continuum limit. 
The present lattice model breaks fourteen of the sixteen supersymmetries at a finite lattice 
spacing. At least, in the perturbation theory, keeping partial supercharges on the lattice 
guarantees the restoration of the full symmetry without the usage of any fine tunings. 
We numerically estimate the supersymmetric Ward-Takahashi identity to verify the SUSY 
restoration beyond the perturbation theory. ' 

3.1 Simulation details 

Numerical simulations are made using the rational Hybrid Monte Carlo method[24]. As 
detailed below, the field configuration of the HMC consists of twelve bosons: two gauge 
fields, eight scalar fields and two auxiliary fields. Moreover, the dynamical effects of the 
fermions are treated using the pseudo-fernrion method with the rational approximation 
with the exception of the phase of the pfaffian. The effects of the phase are included in the 
numerical results of the SUSY WTI presented in this paper, using the phase reweighting 
method. 

The present lattice action has two four-fermi interactions, 

2 2 tr V’-p} 2 ) , (3.1) 

0 t,x p= 0,1 ' ' 

which are of the order of the cut-off. These terms originally come from the two fermi 
terms in the right-hand side of the super transformation of the link fields (2.14). If we 
truncate these terms, the lattice action loses the two exact supercharges that are keys 
in the supersynnnery restoration. Therefore, the terms (3.1) should be included in the 
simulations. By introducing two auxiliary fields, we can express (3.1) as 

54/ = ^- S tr ( a p + VUpkp, V’-p]) • (3.2) 

0 t,x p= 0,1 

Hence, the fermion action is written in a bilinear form which is suited to the simulations. 
Then, to integrate fermions gives the pfaffian, 

pf(-D) = J Df e“ fTOT , (3.3) 

which is generally complex due to the Majorana-Weyl nature of the ten-dimensional fermions 
We deal with the absolute value and the phase of the pfaffian individually. Since det(-D) = 
pf(R) 2 , the absolute value of the pfaffian can be prepared by a complex pseudo fermion (j) 
as 


|pf(£>)| = det(H t H) 1/4 , 


J Dcj) ] D(j) exp|—</) t (R t T))"^} , 


7 See [27] for M = (2, 2) SYM. 


(3.4) 

(3.5) 
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(3.6) 


with the rational approximation of the fourth-root of D^D, 

(Dt D )-V^ ao + g(_^L_). 

To compute the inversions of D'D + fy, we use the multiple shift CG solver [25]. The 
degree N r and the coefficients depend on the accuracy of the approximation and the 

maximal/minimum eigenvalues of D^D [26]. They are determined from the eigenvalues 
at the thermalization step and fixed thereafter. Moreover, we compute the phase of the 
pfaffian for every configurations and use the phase reweighting method to obtain results 
involving the whole dynamical effects of the fermions. 

The classical bosonic action in (2.1) has the following flat directions, 

An = 0, Xj = constant and diagonal matrices. (3-7) 

At high temperature T e g = T /A 1 / 2 2> 1, the HMC stably runs with small fluctuations of AQ 
around zero and Tbi around unity. However, at low temperature T e g < 1, the HMC can get 
unstable by entering a flat direction. Then, the magnitude of the scalar fields R = y] tr(Y 2 ) 
monotonically increases. To avoid the instability, we added a supersymmetry breaking mass 
term[27], 

M _ 9 

Smass = gT- ^ ( 3 - 8 ) 

0 t, x i= 2 

to the action, where mo = ma is the mass in the lattice unit. This controllable supersymme¬ 
try breaking term enable us to investigate the SUSY breaking cut-off effects by examining 
the SUSY WTI and to discuss the restoration of supersymmetry in the continuum limit. 

The trajectory length is 0.5 and the time step of the molecular dynamics is chosen 
to keep the acceptance rate over 80%. We take N r ,oti,/3i of the rational approximation 
to guarantee an accuracy of 10~ 13 and compute the maximum and minimum eigenvalues 
every trajectories for confirmation. 

3.2 SUSY Ward-Takahashi identity 

In the present model, supersymmetry is broken by three sources: finite temperature, the 
mass of the scalars and the lattice spacing. The breaking effects of the temperature and the 
mass are physical, while the effect of the finite lattice spacing is purely a lattice artifact. We 
have to show that the lattice artifact vanishes in the continuum limit. As explained below, 
to observe the symmetry restoration, we can use the Ward-Takahashi identity including 
the breaking effect of the mass term (partially conserved supercurrent) [27]. 

If the action is invariant under supersymmetry, then the SUSY Ward-Takahashi iden¬ 
tity holds as a manifestation of the symmetry. At finite temperature, supersymmetry is 
broken due to the fermions obeying anti-periodic boundary conditions. However we can for¬ 
mally derive the same Ward-Takahashi identity using the path integral formulation. This 
enables us to use the SUSY WTI to examine the behavior of the supersymmetry breaking 
effect which comes from the lattice spacing at finite temperature [27]. 
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In the continuum theory, the SUSY WTI is given by 

9 

777 - 

dpJp(x) = —Y(x), (3.9) 

where fi = 0,1, and the supercurrents are 

N ( i 1 

Jfi = jToTiT/M^oi) + 7^7i7/Ar(TA,W) + (^[X i ,X j ]) [ > . (3.10) 

The supersymmetry breaking mass term gives 

y = N'yMX^), (3.11) 

which is the right-hand side of (3.9). Both the supercurrents and the breaking term Y 
are fermionic operators. 

Let us consider the SUSY WTI as a correlation function because the fermionic one point 
function (J) identically vanishes. We choose, for simplicity, Yp{y) as a source operator. 
Then, the resultant WTI is 

2 

dn(Jn,a(x)Yp(y)) = ^j-( Y a(x)Yp(y)) - <5 2 (x - y)(Q a Yp{y)}, (3.12) 

where Q a are the supercharges. If the SUSY WTI holds, the ratios 

d»(J», a (x)Yp(y)) 

(Ya(x)Y 0 (y)) { > 

should give the dimensionless mass m 2 / A everywhere except at x = y. Therefore, by 
investigating lattice counterparts of the ratios, we can extract the supersymmetry breaking 
effect that comes from the lattice spacing. 

On the lattice, the supercurrent is prepared using only the forward-type operators. 
We use (2.19) for the field strength Toi of (3.10), and (2.12) for the covariant derivative of 
(3.10). Meanwhile, we use the naive symmetric difference operator when calculating the 
gradient of the supercurrent, 

Wi) = £ (3.14) 

M=0d 

3.3 Numerical results 

In this section, we show the numerical results of the SUSY WTI for SU( 2) and at the 
dimensionless temperature T e g = Tj A 1 / 2 = 0.3. At this temperature, for SU( 2), the 
problem of the flat directions, which was explained in section 3.1, occurs. Therefore, 
we added the mass term (3.8) to the action (2.15). Table 1 summarizes the simulation 
parameters. We used four different physical masses and three different lattice spacings. 
The first 6000-10000 trajectories were discarded for thermalization, after that we stored 1 
configuration every 20 trajectories. The other details of the simulations are given in section 
3.1. 
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N t 

N x 

m 2 /X 

Nconf 

12 

6 

1.00 

2000 

12 

6 

0.75 

2000 

12 

6 

0.50 

2000 

12 

6 

0.25 

2000 

16 

8 

1.00 

2000 

16 

8 

0.75 

2000 

16 

8 

0.50 

2000 

16 

8 

0.25 

2000 

20 

10 

1.00 

2000 

20 

10 

0.75 

2000 

20 

10 

0.50 

2000 

20 

10 

0.25 

2000 


Table 1 . Simulation parameters and the number of configurations. The dimensionless temperature 
is T e ff = 0.3 and the lattice spacing is determined by aA 1 / 2 = l/(7V t T e g). 


As explained in section 2.3, the lattice action has many extra vacua. The first inter¬ 
esting result is a confirmation that the presence of the vacua is not an issue in practice. 
For SU( 2), the plaquettes have two minima, +1 and —1. The first one is physical, while 
the second one is a lattice artifact. As can be seen in the histogram of tr(.Poi)/-Y (figure 
1), the plaquettes stay close to one, we never observed any negative values in the present 
parameters set. Moreover, figure 2 shows the contributions of the higher derivative terms 
in the bosonic sector, C (B.16). For all parameters, C is relatively small, only a few percent 
of the action, and approaches zero as a —>• 0. The presence of the mass plays an important 
role in suppressing C because the suppressions become larger as the mass increases. Thus, 
we can conclude that the extra vacua do not affect our main results given below. 

To investigate the restoration of supersymmetry, we compare the lattice counterparts 
of the ratios (3.13) with the expected value m 2 / A. In the present notation of the gamma 
matrices, given in appendix A, the diagonal parts a = f3 show the strongest signal and 
are used for our analysis. We thus have sixteen ratios corresponding to the sixteen su¬ 
percharges, where all the ratios should be equivalent under the 0(8) internal symmetry in 
the continuum limit. The inversion of the Dirac operator is computed for all points-to-all 
points using LAPACK[28]. The correlation functions are averaged over the coordinates to 
reduce the statistical errors. Figure 3 shows an example of the shape of the correlation 
functions. 

Figure 4 shows the ratio for a = j3 = 1 with m 2 /A = 0.25 against t at fixed x = N x /2. 
The errors are estimated using the jackknife analysis. We see that in the middle range of 
the plot, the signal is relatively a constant and in good accordance with m 2 / A. Close to 
the borders (t = 0 ,Nt), we can see deviations from the constant, which correspond to the 
contact term of (3.12). To quantify the recovery of the SUSY WTI, we estimate the value 
of the plateau for each a by performing a constant fit in a rectangular area in which the 
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Figure 1. Histogram of the plaquettes for to 2 /A = 0.25. This includes all plaquettes of every 
configuration. As the lattice spacing approaches zero, the distribution is strongly localized at 

Poi = 1 - 




Figure 2. Higher derivative terms in the bosonic sector C against the lattice spacing. The vertical 
axis denotes C/S' B where S' B is the bosonic part of the action without the auxiliary fields a p . 
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Figure 3. Correlation functions with respect to the WTI. The corresponding parameters are 
N t = 20, to 2 /A = 0.5. All the correlators are presented as a function of t/N t at fixed x = N x /2 for 
a = (3 = 1. The other correlators for a = /3 = 2, • • • ,16 behave in the same way. 


plateau can be clearly seen. The choice of the fit areas is irrelevant to the results as far as 
X 2 /dof ~ 1 . 

Table 2 shows the values averaged over a with x 2 and the fit ranges. We find that the 
fitted values of the plateaus show good agreement with m 2 /A. We also show the results 
as figure 5. Clearly, all the values are on a straight line passing through the origin with a 
slope of one. In the massless limit, our results indicate that the supercurrent is conserved 
as expected. These results tell us that the SUSY WTI holds if we ignore the contact term. 

In the continuum limit, the ratios should be m ?/A everywhere except at the origin. On 
the lattice, they behave as constants only far from the origin because the contact terms of 
(3.12) are smeared. To investigate the behavior of the smeared contact terms, let us count 
the numbers of sites on which the ratios are m ?/A within the errors. To do this, for subsets 
of the lattice points T Q , we define 


where 


x«(r a ) 



(error of R a (x )) 2 ’ 


R / n ^(J M , a (0)Y a (f)) _ nf_ 

a[ ’ (Y a (0)Y a (x)) A ' 


(3.15) 


(3.16) 


We took the largest possible r Q within y 2 (T Q )/A r Q < 1 , where N a are the numbers of 
sites in r Q . Roughly speaking, T a are the plateau areas. In the continuum limit, the 
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Figure 4. as a f unc ti° n of t/N t at fixed x = N x /2 for m 2 /A = 0.25 (denoted by 

the dashed line). 


lattice size 

m 2 / X 

fit value 

X 2 / N do f 

fit range 

12 x 6 

1.00 

1.007 ± 0.043 

0.51 

3x5 

12 x 6 

0.75 

0.750 ± 0.029 

0.59 

3x5 

12 x 6 

0.50 

0.502 ± 0.013 

1.42 

5x5 

12 x 6 

0.25 

0.254 ± 0.006 

1.10 

5x5 

16 x 8 

1.00 

0.990 ± 0.044 

0.59 

5x5 

16 x 8 

0.75 

0.714 ± 0.029 

0.48 

5x5 

16 x 8 

0.50 

0.513 ± 0.015 

0.60 

5x7 

16 x 8 

0.25 

0.249 ± 0.005 

0.93 

7x7 

20 x 10 

1.00 

1.001 ± 0.043 

0.48 

9x5 

20 x 10 

0.75 

0.745 ± 0.026 

0.69 

9x7 

20 x 10 

0.50 

0.482 ± 0.012 

1.03 

11 x 7 

20 x 10 

0.25 

0.247 ± 0.005 

0.93 

11 x 7 


Table 2. Values of the plateau averaged over a. 
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Figure 5. Values of the plateau averaged over a against the SUSY breaking mass to 2 /A. All the 
values are on the dashed line passing through the origin with a slope of one. 



m 2 /?i 

Figure 6. Continuum limit of the values of the plateau obtained for each spinor. The x-axis 
denotes to 2 /A. The red points with errors separated by a little distance denote the sixteen values 
corresponding spinors a = 1, • • • ,16 from left to right. The blue points denote the averaged values. 
All the values are on the dashed line passing through the origin with a slope of one. 
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lattice size 

m 2 / A 

N/V 

12 x 6 

1.00 

0.50 

± 

0.03 

+ 0.07 

- 0.06 

12 x 6 

0.75 

0.49 

± 

0.05 

+ 0.05 

- 0.03 

12 x 6 

0.50 

0.49 

± 

0.06 

+ 0.05 

- 0.06 

12 x 6 

0.25 

0.47 

± 

0.07 

+ 0.05 

- 0.06 

16 x 8 

1.00 

0.57 

± 

0.05 

+ 0.03 

- 0.06 

16 x 8 

0.75 

0.56 

± 

0.07 

+ 0.05 

- 0.04 

16 x 8 

0.50 

0.56 

± 

0.05 

+ 0.05 

- 0.05 

16 x 8 

0.25 

0.53 

± 

0.05 

+ 0.05 

- 0.04 

20 x 10 

1.00 

0.64 

± 

0.04 

+ 0.04 

- 0.05 

20 x 10 

0.75 

0.65 

± 

0.04 

+ 0.03 

- 0.03 

20 x 10 

0.50 

0.65 

± 

0.04 

+ 0.03 

-0.02 

20 x 10 

0.25 

0.63 

± 

0.03 

+ 0.04 

- 0.03 


Table 3. N/V for each parameters set. The first and the second column in N/V denote N/V 
averaged over a and their statistical errors, respectively. While the third column denotes differences 
among spinors. The upper and lower values in the third column represent the maximum and 
minimum values of the differences between each N a /V and the averages, respectively. 


smeared contact terms should become the genuine contact terms. Hence, the ratios N a /V 
must approach one as the lattice spacing decreases. Table 3 shows the N a /V ratios. We 
find that N a /V do indeed increase as the lattice spacing decreases. This means that the 
smeared contact terms tend toward the expected delta function in the continuum limit. 

Finally, we study whether the same results can be obtained for each supercharge. We 
have already seen that the SUSY WTI averaged over a. holds in the continuum limit. All 
we have to do is to see the degeneracy of the results. Figure 6 shows the fitted values of the 
plateau for all the sixteen supercharges. For the extrapolation to the continuum limit, we 
use constant fits. The values coincide with each other within the statistical errors. Namely, 
the ratios are the same as m? /A independently of the spinors. Moreover, N a /V in table 

3 suggest that the contact terms approach the delta function for each spinor. Therefore, 
supersymmetries are recovered for all the spinor indices. 

4 Summary and discussion 

We are performing lattice simulations of two-dimensional SYM with sixteen supercharges 
to verify the duality conjecture. When starting from the lattice action preserving two of 
the sixteen supercharges on the lattice, a non-trivial and crucial question is whether all 
the sixteen supercharges are restored in the continuum limit. The numerical results of 
the SUSY WTI tells us that the SUSY breaking cut-off effect vanishes and the remaining 
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fourteen supercharges are restored in the continuum limit. This is the first result showing 
that the Sugino action does work for two-dimensional Af = (8,8) SYM. 

The lattice Af = (8,8) SYM has many degenerate vacua. However, we have not 
observed the degeneracy in the simulations with the SUSY breaking mass term in SU( 2) 
presented in this paper, nor in subsequent simulations at high temperatures without the 
mass term in SU( 8). When we start the simulations from the trivial configuration (U^ = 
1 ,Xi = 0), at high temperature T/ A 1//2 2> 1, both temporal and spatial Wilson loops and 
the plaquettes tend to stay close to unity, and the scalar fields distribute around zero. In 
other words, transitions from trivial vacuum to the others would be suppressed at high 
temperature, for sufficiently large N, or near the continuum limit so that our choice of 
the lattice held tensor (2.19) is justified. The transitions were, somewhat surprisingly, not 
observed at low temperature such as T/X = 0.3 in IV = 2 and at the lattice spacing 
aX 1 / 2 ~ 0.28 at which the SUSY WTI were measured. The degenerate vacua might not be 
a problem for a wide range of simulation parameters. 

The next step of this study is to investigate the phase structure and to test the duality 
conjecture[18]. To see the correspondence between the gauge theory and the black 1- 
branes, we have to choose simulation parameters in which the theory is well-described by 
the black 1-branes. Then, from lattice simulations with the determined parameters, we 
will estimate the internal energy of the branes[22]. To do such verification of the duality, 
further simulations are ongoing. 
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A Definition of twisted fields 

We use two different notations of the fields in this paper. The action (2.1) is written 
in the original variables for /i = 0,1, i = 2, ••• ,9, a = 1, • • • , 16, while the 

twisted action (2.5) is written in the twisted variables, A^, Bi,C,<p± and r /±, . x±i for 

H = 0,1,2,3 and i = 0,1,2. In this appendix, we give the relation between the two 
notations. 

Let us express the action (2.1) written in as the Q±-exact action (2.5) 

written in the twisted variables A^, Bi, C, <f>±. The gauge fields are unchanged, and the 
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other bosonic twisted fields A 2 , A3, -Bj, C, 4 >± are given as follows: 


— X (/i — 2,3), 

Bo = Xq 1 
Bi = X S , 

B 2 = -X 4 , (A.l) 

C = 2 X 7 , 

<t>+ = Xg — *Ag, 

4>- = Xg — i-X’g. 

The fermionic twisted variables 77 ±,V , ±/z>X±i and can be related to each other by 


/ T 1 \ 


1 V’-O + ^ 

T 2 


V'-l - *X +2 

^3 


i’-2 + *x+i 

t 4 


V ’-3 - *X+0 

^5 


V’+o + 5 V- 



V'+l + * X -2 

t 7 


-0+2 - ix-i 


1 

" 0+3 + * X -0 

Tg 

“71 

1 

1 

0 

1 

bO|^. 

C3 

+ 

®10 


+ ix+ 2) 

T11 


-ity -2 - *X+i) 

®12 


-*(^-3 + * X + o ) 

®13 


1 

(S>. 

+ 

0 

1 

LO|^- 

1 

T 14 


-*(</>+1 - *X—2) 

^15 


- i (^+2 + *X-i) 

\® 16 / 


^-*(^+3 - * X - o )/ 


with the following representation of gamma matrices, 

71 = cr 2 <g> 1 (8) cr 3 (8) (72, 

72 = 02 ® 1 < 8 > cr 2 < 8 > 1 , 

73 = 02 <g) 1 <g) (71 <g) ( 72 , 

74 = (J 3 (g) ( 7 2 (8) 1 <g> (72, 

75 = 03 <8> (7 2 <g> <7 2 <8> <73, (A.3) 

76 = <7 3 (g> <7 2 <g> <7 2 <g> <7i , 

77 = <73 (g> <7i <g> 1 <g) 1 , 

78 = (73 <8) <73 ® 1 (g) 1 , 

79 = (71 <g) 1 <g> 1 <g) 1, 
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where 07 are the Pauli matrices. 8 

The two scalar supercharges Q± are defined as 


Q+ 2Q13), 

Q- = ^(Qi _ iQv)- 


(A.5) 

(A. 6 ) 


The Q -(--transformations of the fields come from those of Q a given in ( 2 . 2 ), (2.3) and (2.4). 

Using the definitions of the twisted fields above and the (5±-transformations (2.8), and 
adding the Gaussian integrals of the auxiliary fields and Hi to the action (2.1), we can 
rewrite the action ( 2 . 1 ) as the Q±-ex&ct form (2.5). 


B Lattice action written in twisted variables 

Here we give the full lattice action by performing the (^-transformations (2.14) in the 
lattice Q±-ex act action (2.15). First, let us introduce some useful notations and covariant 
difference operators. Then, by using them, we present the lattice boson and fermion actions. 


The following shifts 

tp +l *(x) = U ll (x)(p(x + p,)Ul(x), (fj, = 0,1), (B.l) 

= u l(x- p,)<p(x- p)U,j(x- fi), (M = 0,1), (B.2) 

p ±fl (x) = <p(x), (n = 2,3) (B.3) 

and covariant difference operators 

V±<p = ±{<p k »-<p), (/z = 0,1), (B.4) 

Vj v = i[A li ,v\, (n = 2,3) (B.5) 

are useful to write the full lattice action. Here, the sites on which the fields are defined 
were abbreviated. This is possible because the gauge covariance always tells us where the 
fields live. The lattice action has two complex covariant difference operators, 

V+*V = + P^ +tl - VP*? ~ P^), (B.6) 

= \{P V & + <pP„v - P^^ ~ V-"P^\ (B-7) 

for /i, v = 0 , 1 , n 7 ^ v. 


8 The gamma matrices act on the fermions ip a as ( 7 i) aiQ2 'I' a2 where the spinor index a corresponds to 
the indices of the Pauli matrices as follows: 

(7l)aia 2 = (<72)ui2 ® lji j 2 ® ® (^hh, (A A) 

with a — 8(j — 1) + 4 (j — 1) + 2 (k — 1) + l. 
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After performing the lattice (^-transformations (2.14), we obtain the full lattice boson 
action, 

Sb = ^^trj^ + ,(/>_] 2 + \[C,4>+][CAA ~\[C,Bi} 2 

t.X ^ 


-V+<£ + V+^_ + [Bi^+WBi,#-] + ^(V+G ) 2 (B. 8 ) 

+{Hi + itpi) 2 + $ + + iG^) 2 + G 2 |, 

where pi and G M are given by 

<p 0 = V+A 3 +ViA 2 -i[B l ,B 2 ], (B.9) 

<Pi=ViAa-VoA 2 -i[B 2 ,Bo], (B.10) 

ip 2 = i[A 2 ,A 3 ] + F m — i[B 0 , Bi], (B.ll) 

and 

G 0 = i[A+°, B 0 ] - i[A 2 , B+°] + VfB 2 , (B. 12 ) 

G, = i[A+\ B 0 ] + i[A 3 , B+ 1 ] - Vq 1 S 2 , (B.13) 

G 2 = -V^B 0 + V+B l + i[A 3 ,B 2 ], (B.14) 

G 3 = — V 0 Bo — VfBi — i[A 2 , B 2 \. (B.15) 


The lattice boson action is semi positive definite. The actions of the seven auxiliary 
fields and Hi are given by the Gaussian integrals which can be analytically integrated. 
Every terms except for ipi and G^ look like the continuum ones. For pi and G^, some fields 
are shifted. Since the higher derivative terms in the bosonic sector, 

C =y o Yl tr {*^01 [A 2 , A 3 ] - Vq A 2 VrA 3 + V+A 3 V+A 2 
t,x —iFoi[Bo, Bi] + V 0 B^V^Bi — VqBiV 1 Bq 

—iV 0 A 2 [Bq, B 2 ] + zV 0 Bq[A 2 , B 2 ] — zV 0 1 B^A^ 1 , Bo] 

-iV+A 3 [B 1 ,B 2 ] + iV+B 1 [A 3 ,B 2 ]-iVo 1 B 2 [A 3 ,B+ 1 ] (B.16) 

-iWtA 2 [B 1 ,B 2 ] + iV+B 1 [A 2 ,B 2 ]-iVfB 2 [A 2 ,B+°] 

+iV^A 3 [B 0 , B 2 ] - iV±B 0 [A 3l B 2 ] + iVfB 2 [A+°, B 0 ] 

+ [A^A^B,] - [A+ 1 ,5o][^3,B 1 +1 ] + [A 2 ,B+°][A+°,Bo]} , 

disappear in the naive continuum limit, the lattice boson action reproduces the continuum 
boson action in (2.1) by using (A.l). 
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The fermion action is given by 

Sf = 2 "^^ tr { +X+il<t>-,X+i] +^+m[s(^- + 0 ^)>' 0 +m] 

+ ^-[0+,??-] +X-»[0+.X-i] +^-#i[i(^+ + 0+ i ),'0-p] 

- ^+[C,r/_] +x+i[C,X-i] + ^+AU C + C+fi )^-^ 

~ V-[ B i,X+i] ~ r]+[Bi,X-i\ + 2 e ijk X-i[Bj,x+k] 

~ iW ^rj+t/)-^ ~ iVnV-ip+iJi + W-, V+W-n + [A M , V-W+n 

- 2 x+o(+*V ( j'' 0 _ 3 + iVfip - 2 + [A^" 1 , V’-i] + [^3"°> V’-o]) 

+ 2x-o(+*V(}~V ; + 3 + iV ] ^'0_|_2 + [A^ 1 , V’+i] + [Aij~°, V’+o]) 

- 2 x+i(-*Vq ^-2 + iVf^_ 3 - [Aj 0 ,^] + [A^VCJ]) 

+ 2 x-i(-*Vq ^ +2 + *Vpi /> +3 - [A 2 0 , V^+o] + [A 3 

- 2x+2(+*V(j'V~i - ?'V|V-o - [A 2 ,^-3] + [-4.3, -0—2]) 

+ 2 x-2 (+*V([ 1 '0 +i - *V+>+o - [A 2 ,'0+3] + [-4.3, -0+ 2 ]) 

- 2 ^- 3 [V>Vvo°] - 2^-2 [^o -1 ^;! 1 ] 

- 2 ^- 3 [ j B 1 +1 ,^+i] + 2 ^_ 2 [ J B+°, V’+o] 

+ 2^ +3 [B 0 -°, VC 0 °] + 2^+2 [B^ 1 , VC1 1 ] 

+ 2 ^ +3 [ 5 i +1 ,^-i] - 2 ^ +2 [-B+ ( \V’-o] 

- 2^-3 [B 2 , V^+2] + 2^+3 [B 2 , V>-2] 

- #-o[A^"°, [Bo,i) + o]\ - #_o[So, [A^^+o]] 

- iip-o[[A 2 , [Bf°,i> + 0]] - [A 2 ,^+o]] 

+ iV»-1 [A3, [-B+ 1 , V’+i]] + [A 3 , i/} +1 ]] 

- iV’-iIAj' 1 , [-Bo, V’+i]] - #-i[B 0 , [- 4 + 1 , - 0 +i]] 

+ Cp 

- \ H w+p^-p} 2 }, 

p=o,i J 

re 

£p i ^ 10 ^ 2 ~ -B2-P01)” 1 + -B2-P10 — -P01-B2, V’+o} 

+ V>-o-Poi^o s 2 - ^-o-B 2 Y’+o p io + ^lo5 2 -i/> + o-Poi - ^o-Pio^+o-B 2 

- 2^” 1 {(-B0I-B2 - ^2-Plo) -1 + B 2 Poi — P10B2, V , +l} 

- V>-i-Pio^?-B 2 + ^-i-B 2 V>|?-Poi - ^^2^+iPio + 1 >±%PoiiJ>+iB 2 
+ ^-0^01^+1-82 - V’-o^V'+i-Pio + ^-o-Blo^^+i - ^- 1 ^+iB 2 Poi 

- i>-o^ti p oiB2 + ^-0-B 2 PlO^+? + ^B 2 ^Pqi - ^tlPl^Xl B 2 

- ^_iPio^+o-B 2 + ^-i-B 2 V’+o-Poi - VAi-Poi-^Y’+o + ^V’+o^-Pio 
+ ^-i^o-Pio-Bs - V’-i-BaPoi^d - ^-iB 2 ^toPio + ^P m ^Xl B 2 - 


(B.17) 


(B.18) 
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In the naive continuum limit, the complicated term Cp becomes two simple terms, 

C P = -2^-i [B 2 ,ip+o] + 2^- 0 [-B 2 ,^+i], (B.19) 

while the shifts of fields in (B.17) disappears. Thus, by using (A.l), (A.2) and (A.3), we 
find that the lattice fermion action reproduces the continuum fermion action in (2.1). 
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